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We report experimental determination of the phase-matching function for type-II three-wave 
mixing in a periodically poled KTiOP04 waveguide in the 792-815 nm spectral region. The mea- 
surement was performed by sum- frequency generation of spectrally tuned fundamental components. 
Strong dependence of the observed signal on the excited spatial modes in the waveguide has been ob- 
served and fully interpreted. These results indicate a route to employ the waveguide for spontaneous 
parametric down- conversion producing photon pairs in well-defined spatial modes. 
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Photon pairs generated in spontaneous parametric 
down-conversion find multiple applications in the emerg- 
ing field of quantum-enhanced technologies 1 3 . Standard 
sources of photon pairs based on bulk crystals suf- 
fer from low brightness and complicated modal struc- 
ture, which hinders their utilization in practical systems. 
These features can be improved substantially by realiz- 
ing the down-conversion process in nonlinear waveguides. 
The spatial structure of the generated photons is then 
defined by the waveguide, and the confinement of the in- 
teracting fields to a small transverse area throughout the 
waveguide length enhances the efficiency of the down- 
conversion process. In addition, waveguide manufactur- 
ing can include the technique of quasi-phase matching 
(QPM), which allows one to employ larger nonlinear co- 
efficients that are otherwise inaccessible at desired wave- 
lengths. All these benefits have been demonstrated in 
recent experiments using periodically poled potassium ti- 
tanyl phosphate (PP-KTP) waveguides 4 7 . 

Typically available PP-KTP waveguides support sev- 
eral transverse modes in the 800 nm wavelength region, 
where photon pairs have been generated experimentally. 
This may have a deleterious effect on the modal purity 
of the generated photons, if the down- conversion process 
takes place simultaneously to different combinations of 
spatial modes. However, the phase matching condition 
in a waveguide depends on propagation constants, which 
differ between the modes due to intermodal dispersion 8 . 
This dependence offers a chance to single out a partic- 
ular down-conversion process into a well-defined pair of 
transverse modes, thus delivering spatially pure photons. 

In this Letter we present an experimental determina- 
tion of the phase matching characteristics for a type- 
II three-wave mixing process in a multimode PP-KTP 
waveguide. This provides detailed spectral and spatial 
information that is instrumental in designing the down- 
conversion process with desired modal properties. The 
measurement was carried out using the inverse process 
of sum- frequency (SF) generation with spectrally tuned 
fundamental beams. 

The experimental setup is presented in Fig. 1. The 
sample was a 4.8 mm long PP-KTP structure (AdvR 
Inc.), which contained a number of ion-exchanged waveg- 
uides with a varying poling period, suitable for type- 



II three-wave mixing with fundamental beams in the 
800 nm wavelength region. A selected waveguide was 
excited with a beam that combined two orthogonally po- 
larized components. The light source was a home-made 
Ti: Sapphire femtosecond oscillator (800 nm central wave- 
length, 25 nm FWHM bandwidth) followed by an optical 
isolator OI. The laser light with polarization rotated to 
45° with respect to the plane of the setup was sent to a 
Michelson-type interferometer. The input beam was split 
using a polarizing beam splitter PBS1 into two compo- 
nents, whose spectra were tuned individually with nar- 
rowband (approx. 0.6 nm FWHM) interference filters IF 
mounted on rotation stages. The spectral range between 
792 nm and 815 nm was covered with two sets of filters, 
designed for 800 nm and 815 nm. The selected filter 
bandwidth offered a good compromise between the spec- 
tral resolution of the measurement and the level of the 
sum-frequency signal. In order to extract the combined 
beams from the interferometer their polarizations were 
rotated by 90° with quarter- wave plates A/4 placed in 
both arms. The piezo-driven mounts for end mirrors PM 
were used to correct for the varying deviation of the out- 
put beam introduced by rotating the interference niters. 
The feedback signal to control the tilt of the mirrors was 
obtained from a 4% pick-off beam splitter BS1 that sent 
a beam to two quadrant detectors QD. The pick-off beam 
was also used to monitor the wavelength of the interfer- 
ometer output with a spectrometer SP. 

The prepared beam was focused at the input face of 
the waveguide using an aspheric lens AL with 8 mm fo- 
cal length. The output light was collected by a birefrin- 
gence free, infinity corrected 50 x microscope objective. 
Dichroic mirrors DM and color filters F were used to 
separate the SF signal from the fundamental. The in- 
tensities of the two fundamental components, separated 
on a polarizing beam splitter PBS3, were measured us- 
ing silicon photodiodes PD1 and PD2, while the SF sig- 
nal was directed to a GaP photodiode PD3 connected 
to a lock-in amplifier. It was ensured that the detectors 
collected only light transmitted through the waveguide. 
Flipper mirrors were used to alternatively direct the out- 
put beams to CCD cameras in order to record magnified 
images of transverse intensity distributions at the exit 
surface of the waveguide. 
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FIG. 1: Experimental setup. See text for symbol explanation. 




FIG. 2: Transverse intensity distributions of exemplary modes 
observed in the waveguide at (a) fundamental and (b) sum- 
frequency wavelengths, obtained by imaging the output plane 
of the waveguide onto a CCD camera. 
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The advantage of the constructed setup was indepen- 
dent control of input coupling for both the components 
of the fundamental beam by tilting the end mirrors PM. 
This enabled us to selectively excite specific transverse 
spatial modes, which was verified with images collected 
on the CCD camera, as depicted in Fig. 2(a). We la- 
bel the waveguide modes with two integers ij specifying 
the number of nodes in directions parallel and perpen- 
dicular to the top crystal surface and an index V or H 
corresponding to the vertical or horizontal polarization of 
the fundamental component, or S for the sum- frequency. 
We observed that the SF signal depended critically on 
the coupling of the fundamental components into specific 
spatial modes. The SF signal was typically generated in 
a well defined spatial mode, as shown in Fig. 2(b). These 
features indicate the important role of intermodal disper- 
sion in the phase matching of three- wave mixing process. 

More detailed information about the phase matching 
can be obtained from a scan of the SF signal intensity as 
a function of the wavelengths of the fundamental com- 
ponents. In order to aid interpretation of the data, we 
aimed at exciting well defined spatial modes of the fun- 
damental components throughout the scan. For this pur- 
pose, before each scan the coupling of the fundamental 
components was adjusted with the end mirrors PM to 
optimize the intensity of the SF signal at a chosen wave- 
length. The coupling was then kept fixed during the scan 
by locking the feedback loop of the piezo drivers to the 
selected position. The measured sum-frequency intensity 
was normalized by dividing it by the product of the in- 
tensities of the fundamental components measured after 
the waveguide. This compensated for variable transmis- 
sion of the interference filters and enabled comparison 
between scans obtained with different coupling efficien- 
cies. 

The measured maps of the normalized SF intensity are 
shown in Fig. 3. It is seen that maps are composed of a se- 



FIG. 3: The measured normalized sum- frequency intensity 
vs. wavelengths of horizontally and vertically polarized fun- 
damental components (a,b), scanned with 0.2 nm step, and 
corresponding one-dimensional cross sections at degenerate 
wavelengths (c,d). The components were optimized for cou- 
pling in (a,c) fundamental 00 modes and (b,d) a combina- 
tion of 00 v and 02#. Each band is labeled with a triplet 
ijv + i' Jh kls of interacting spatial modes, identified fol- 
lowing the method described in the text. Imperfect excita- 
tion of selected spatial modes is seen to generate additional 
bands. The dashed lines in (a) indicate the spectral region of 
the down- conversion process defined by a pump field centered 
at 399.8 nm with 1 nm FWHM. 



ries of bands. Each band can be associated with a triplet 
of spatial modes — two fundamental with orthogonal po- 
larizations and one SF — for which the phase matching 
condition is satisfied. The identification of the relevant 
modes turned out to be non-trivial due to spectral over- 
lap of the bands and impure excitation of the spatial 
modes. The starting point of the identification proce- 
dure was a series of one-dimensional scans with identical 
wavelengths for both the fundamental components, sim- 
ilar to those shown in Fig. 3(c,d), carried out for various 
spatial couplings. This provided unambiguous informa- 
tion on bands that were the brightest or well isolated in 
the spectral domain. 

Further identification required a theoretical model. 
The phase matching condition for a specific triplet of 
modes ijv + i' j'jj kls in the waveguide is given by the 
equation: 

2ttc 

uvn¥j(uv)+UHnvj'(uii) = (wy+cj#)ra£ z (cjy+cjtf )- — 

(1) 

where ujy and ujh are the frequencies of the two funda- 
mental components, and n(-^ri^-, , and are frequency- 
dependent effective refractive indices of the relevant 
modes. The parameter A is the effective grating pe- 



3 



8 5 T i 8 t 



**, , ^ up, ** ^W*,v ay, ,?fl; , Hi , 

794 796 798 800 802 804 806 808 810 812 814 

Wavelength (nm) 



FIG. 4: Triplets of the interacting spatial modes ordered according to the wavelength of a frequency-degenerate process. Each 
triplet is represented by a pair of points indicating the measured (+) and calculated (x) sum- frequency intensity, normalized 
to 100 for the OOy + 00h — > 00s band. Triplets involving 00s mode are shown in black, while those involving higher SF modes 
in red. 



riod of the QPM structure. Refractive indices can be 
decomposed 8 into the material and geometric contribu- 
tions according to nfj(uj) = n v (uu) + AnJ^, and similarly 
for H and S. Here n v \uj) is the refractive index of bulk 
KTP and the geometric correction Anfj is assumed to be 
constant within the frequency range swept in the experi- 
ment. Given two identified triplets that differ only by one 
of the participating modes, for example ijy -hi' ' j' H kls 
and ijy + i' j' H «->• k'l' s , this assumption allows us to de- 
termine the difference An^ — Anf, z , . Using bulk KTP re- 
fractive indices 9 and data obtained from one-dimensional 
scans, we were thus able to calculate relative spectral lo- 
cation of other pairs of bands that involved Anf z and 
Anf,j,. These predictions matched closely the locations 
of bands revealed by one-dimensional scans, with devia- 
tions less than 0.2 nm. The procedure yielded a complete 
identification of triplets, as indicated in Fig. 3(a,b). 

The intensity of the sum-frequency signal is propor- 
tional to the squared spatial overlap of the triplet of 
the modes involved 6 . We calculated electric field dis- 
tributions of the modes using a vector finite- difference 
modesolver 10 , assuming error function depth profile with 
refractive index modulation taken from Roelofs et al 8 
The geometric parameters of the waveguide were ob- 
tained by coupling in white light and analyzing the im- 
age of the output face of the waveguide. The theoreti- 
cally calculated squared overlaps are compared in Fig. 4 
with normalized intensities of the sum-frequency signal 
measured for degenerate wavelengths of V and H com- 
ponents. In this measurement, for each selected band 
the intensity was maximized by manually adjusting the 
input coupling. A reasonable agreement between mea- 
sured and simulated data is seen. Discrepancies can be 
attributed to impure mode excitation and imperfect nu- 



merical modeling. 

The experimental data presented in this Letter pro- 
vide detailed information on the role of transverse spa- 
tial modes in three-wave mixing, which can be used to 
engineer the process of spontaneous parametric down- 
conversion. Suppose that the pump field S is launched 
into the waveguide in the fundamental 00s mode. Using 
Fig. 4 it is easy to identify the pairs of modes in which 
down-converted photons can be produced, and the lo- 
cation of the respective spectral bands in the (Ay,A#) 
plane. Thanks to intermodal dispersion, the center of the 
00y + 00 h 00s band is shifted from other combina- 
tions by more than 3 nm. Therefore the 00s —> 00y +00# 
process can be separated by using a long waveguide, for 
which the bands become sufficiently narrow. Further, the 
wavelengths of the generated photons must lie within a 
region defined by the bandwidth of the pump field, as 
marked in Fig. 3(a) with dashed lines. Because the group 
velocities for V and H polarized photons are unequal, 
the slope of the phase matching bands will differ from 
the slope of the region defined by the pump field. As 
a result, photons generated in different pairs of spatial 
modes will be separated spectrally. This should allow 
one to isolate 00y + 00h pairs by employing spectral fil- 
tering only. Note that the spatial purity of the pump field 
mode is critical in this scheme, as in the vicinity of the 
00s ^ 00y + 00 h band there are other bands involving 
higher-order S modes. Thus we see that the interplay be- 
tween spatial and spectral degrees of freedom provides a 
wealth of tools to engineer the down-conversion process. 

We acknowledge helpful discussions with LA. Walm- 
sley and C. Silberhorn. This work was supported by 
EU FP6 project QAP (Contract no. 015848) and MNiSW 
grant no. N N202 1489 33. 



* Electronic address: mkarp@fuw.edu.pl 

1 A. V. Sergienko and G. S. Jaeger, Contemp. Phys. 44, 341 
(2003). 

2 I. Marcikic, A. Lamas- Linares, and C. Kurtsiefer, Appl. 
Phys. Lett. 89, 101122 (2006). 

3 P. Trojek and H. Weinfurter, Appl. Phys. Lett. 92, 211103 



(2008). 

4 K. Banaszek, A. B. U'Ren, and I. A. Walmsley, Opt. Lett. 
26, 1367 (2001). 

5 A. B. U'Ren, C. Silberhorn, K. Banaszek, and I. A. Walm- 
sley, Phys. Rev. Lett. 93, 093601 (2004). 

6 M. Fiorentino, S. M. Spillane, R. G. Beausoleil, T. D. 



4 



Roberts, P. Battle, and M. W. Munro, Opt. Express 15, 
7479 (2007). 

7 A. Trifonov and A. Zavriyev, J. Opt. B 7, 772 (2005). 

8 M. G. Roelofs, A. Suna, W. Bindloss, and J. D. Bierlein, 
J. Appl. Phys. 76, 4999 (1994). 



9 K. Kato and E. Takaoka, Appl. Opt. 41, 5040 (2002). 
A. B. Fallahkhair, K. S. Li, and T. E. Murphy, J. Light- 
wave Technol. 26, 1423 (2008). 



